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Q_i, Abstract 

A systematic study of optical and transport properties of the Hubbard model, based on 
Metzner and Vollhardt's dynamical mean-field approximation, is reviewed. This model 
shows interesting anomalous properties that are, in our opinion, ubiquitous to single- 
band strongly correlated systems (for all spatial dimensions greater than one), and also 
compare qualitatively with many anomalous transport features of the high-T c cuprates. 
This anomalous behavior of the normal-state properties is traced to a "collective single- 
band Kondo effect," in which a quasiparticle resonance forms at the Fermi level as the 
temperature is lowered, ultimately yielding a strongly renormalized Fermi liquid at zero 
temperature. 
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1 Introduction and survey 





Figure 1: Typical experimental results for resistivity (p), NMR relaxation rate (T\), Hall coef- 
ficient (Rh), and inverse mobility (l/n) in high-T c compounds (from j^j). 

The discovery of the high-T c superconductors based on CuO-compounds has led to a 
large amount of theoretical work about the peculiar properties of these materials. While initial 
work concentrated on explaining "why T c is so low" and possible exotic mechanisms for the 
superconductivity, it became obvious that an understanding of the superconducting mechanism 
is linked to an understanding of the anomalous normal-state properties of these compounds || . 
Most prominent among these are the linear (in T) resistivity, a linear (in T) NMR-relaxation 
rate of the Cu-spins, and a Hall angle that grows like T 2 over a rather wide temperature 
region (see Fig. [j] ||). Furthermore, the optical conductivity shows a Drude peak with a width 
1/r ~ T consistent with the linear behavior of the resistivity and a pronounced mid-IR 
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Figure 2: Experimental optical conductivity of TlBa2Ca2Cu^Og (from j^J). Note the prominent 
mid-IR peak at about 1500 cmr x . 



peak at frequencies above the Drude peak (see Fig. |] 0). 

Early on, it was argued || that most of these anomalous properties can be explained by two 
special features appearing simultaneously in these materials: (i) They are strongly correlated, 
i.e. their (effective) local Coulomb interaction is comparable to or larger than the characteristic 
kinetic energy of the relevant carriers, and (ii) they are highly anisotropic with the electrons 
confined to the CuO-planes characteristic for these compounds. This mixture of strong correla- 
tions and effective 2D-character is unusual and can only be found (apart from the cuprates and 
other perovskites) in some Uranium-based heavy-fermion compounds. The great interest in the 
high-temperature superconductors has led to a number of interesting new theoretical ansatzes, 
that, although based on the assumption of strongly correlated carriers, focussed mainly on the 
2D-character of the CuO-planes. Most prominent are Anderson's concept of resonating valence- 
bonds |7| , his recent efforts to understand the electronic properties phenomenologically in terms 
of a Luttinger-liquid fixed-point (known from the theory of ID electron systems) [§, || , Varma's 
phenomenological marginal Fermi-liquid theory [|Hj], and Lee's anyon concept emerging from a 
mapping of the 2D electron system onto a nonlinear sigma- model [|IT|. While the second and 
third ansatz have not been ruled out, both the first and the last are insufficient to account for 
the experimental data of these systems. 

However, even in Anderson's and Varma's theories, one needs additional assumptions about 
relaxation times [[| or dynamical quantities [|10j to consistently account for all anomalies. Both 
theories also lack a thorough microscopic foundation. It is unclear how the special form of the 
single-particle self energy and the dynamical susceptibility (necessary for the marginal Fermi- 
liquid picture) will emerge. Even introducing a special three-body scattering vertex |T^] does 
not resolve this difficulty, since it merely shifts the problem to another level. A foundation for 
Anderson's Luttinger-liquid picture is even less firm. The special feature of ID systems — that 
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the exchange of charge over a given site is strongly hindered by the Coulomb correlations, while 
the exchange of spin costs almost no energy |]TB[ (which leads to the strict separation of spin and 
charge degrees of freedom) — may already be absent in 2D, since a physical electron may move 
"around" any other via a path in the plane. Recent studies on the stability of the possible low- 
energy fixed points in interacting fermion models also suggest that it is the Fermi-liquid, rather 
than the Luttinger-liquid fixed point, that is stable in 2D (unless the small-momentum-transfer 
couplings become singular) |T4[ . 

There is another feature of strongly correlated systems that is usually ignored in dealing with 
the cuprates. From the experimental and theoretical investigation of heavy-fermion compounds 
it is known that the strong correlations themselves already lead to anomalous features in the 
magnetic and transport properties of metals. These systems are usually three-dimensional 
in nature and, despite their complicated bandstructure one can explain their physical 
properties qualitatively (and to some extend even quantitatively) by merely taking into account 



one strongly correlated, localized band hybridizing with one uncorrelated conduction band [16 



This suggests that the prominent anomalous features of the heavy-fermion materials are intrinsic 
to the strong local correlations. Thus to unambiguously understand the origin of any anomalous 
behavior in systems with strong electronic correlations, it is necessary to discriminate between 
those properties that are intrinsic and arise directly from the strong (local) electron-electron 
scattering, and those that are connected to the geometry of the underlying lattices, e.g. the 
2D-character of the Cu-0 planes in the cuprates. 



In previous publications [117L 18 1 we demonstrated that even for the simplest model of strong 



local correlations, the single-band Hubbard model, interesting anomalous results for several 
different transport quantities occur which appear to be intrinsic properties of strongly correlated 
systems. This latter statement is motivated by the fact that our calculations use Metzner 
and Vollhardt's "dynamical mean- field approximation" which may be viewed as the Fermionic 
equivalent of the standard non-local mean-field description of models for magnetism. As is well 
known, such a "proper" mean-field description will in general reproduce the intrinsic features 
of the underlying model. Fluctuation corrections should only affect phase transitions and the 
quantitative values of temperature scales. 

In this review we present a more detailed discussion of the transport properties of strongly 
correlated systems described by the Hubbard model in Eq. (|T|). Special emphasis is made to 
thoroughly discuss the variation of quantities with temperature and doping. We also try to 
make qualitative contact with recent experiments on optical and transport properties of the 
cuprates. The paper is organized as follows: In the next section an introduction into the 
underlying theory (the so-called dynamical mean-field theory) is given. We also give a brief 
account on how the transport properties are calculated. In section 3 the systematics for several 
quantities as functions of temperature and doping are presented and compared to experiment. 
A discussion and outlook in section 4 conclude the review. 
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2 Theoretical background 



The model discussed in this contribution is the single-band Hubbard model [19 



he 

ka 



Our notation is the following: c- is the creation operator for an electron in a Bloch state 

ka 



with wavevector k and z-component of spin a, and n ia is the electron number operator for 
electrons localized at lattice site i with spin a (i.e., the number operator for electrons in Wannier 
orbitals centered at lattice site %). For the sake of simplicity we restrict ourselves to a simple 
hypercubic lattice (in <i-dimensions) with lattice constant a and nearest-neighbor transfer t, 
i.e. = — 2tY,i=i cos(kia); U denotes the local Coulomb repulsion. The energy scale is set by 
Adt 2 =: t* 2 and the choice t* = 1 (this is the required scaling for nontrivial results as d — > oo). 

The Hubbard model is perhaps the simplest model of a correlated electronic system, since 
it consists of a single band of delocalized electrons subject to a local Coulomb interaction. It is 
easier, therefore, to discriminate intrinsic effects of the correlations (induced by the Coulomb 
interaction) from bandstructure effects. Obviously, a more realistic description of a condensed- 
matter system must take into account the existence of more than one band, hybridizations 
between different bands, and also additional local and long-range Coulomb interactions. How- 
ever, these more realistic models may be mapped onto the Hubbard model (|l|) when excitations 



into these other bands occur only virtually, leading to renormalized values for e.g. p0| . Since 
the purpose of the present contribution is to identify effects intrinsic to the local correlations, 
these other interactions may be neglected. 

Models of interacting electrons on a lattice are difficult to solve because they include both 
local and nonlocal parts in the Hamiltonian. The fundamental quantum-mechanical princi- 
ple of complementarity implies that one cannot expect to describe such a system by either 
purely localized or purely delocalized states. Approximations starting from either side (i.e. 
bandwidth/?/ C 1 or bandwidth/?/ ^> 1) are not generally applicable to the intermediate 
regime (bandwidth « U). To obtain sensible results for the Hubbard model in this intermedi- 
ate regime, the Hamiltonian (|l|) is usually studied with exact diagonalization techniques [21 



or quantum Monte Carlo methods (QMC) P2" |. Both methods share the problem that they 



can only be applied to comparatively small system sizes, (~ 20 sites in exact diagonalization 
and ~ 100 sites in QMC [for not too low temperatures]). This restriction seriously affects 
any attempt to determine low-energy/low-temperature properties and makes the calculation of 
dynamical quantities and transport properties difficult. An ansatz is therefore needed to calcu- 
late different physical quantities of strongly correlated systems that (i) maintains the important 
local correlations and (ii) allows calculations to be performed in the thermodynamic limit. 
Such a method was proposed by Metzner and Vollhardt |23| and Miiller-Hartmann f21 | 



who observed that the renormalizations due to local two-particle interactions like the Hubbard- 
U become purely local as the coordination number of the lattice increases. More precisely, 
the irreducible single-particle self energy ^t(z) and the irreducible two-particle self energy 
^k k+qy z i z ') both become independent of momentum for large coordination number (2d — > oo) 
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S|2J: 

Jim E n (z) = E(z) , 

(2) 

Jim r £ g +ir (z, z ) = r(z, z ) . 

One can use standard techniques of field theory to show that the second relation in Eq. (0) 



is a necessary consequence of the first for all two-particle self energies |26|]. A further 
consequence of the first relation is that the solution of the Hubbard model may be mapped 
onto the solution of a local correlated system coupled to an effective bath that is self-consistently 
determined by the following five-step process |27|, |28], |2§], |3U], [JTJ : 



• (i) Choose a suitable starting guess for the local self energy E(z). 

• (ii) Calculate the local single-particle Green function from its Fourier transform 

G-^z) = 1 V— , (3) 

with Gr := l/{z + fi — t£) the noninteracting momentum-dependent Green function. 

• (iii) Obtain the effective medium, denoted by Q(z), by subtracting off the local correlations 
from the Coulomb interaction at site i, 

[G{z)]- 1 :={G^- l + ^) • (4) 

• (iv) Solve the local impurity problem defined by the effective medium Q{z) and the local 
Coulomb interaction Un^riii a ^ s ^ e * ( v ^ a an exac ^ quantum Monte Carlo algorithm, or 
an approximate algorithm based upon a diagrammatic analysis) to obtain a new local 
single-particle Green function G(z). 

• (v) Obtain a new self energy via E(z) = — [^(-z)] -1 and repeat steps (ii) through 
(v) until the local Green functions are identical between two successive iterations [Ga(z) = 
G(z)}. 

This mapping onto a local problem coupled to an effective bath is reminiscent of the standard 
mean-field description of spin systems. As is well known, such a theory emerges systematically 



from the same limit of large coordination number fl32| . One major difference between electronic 
models and spin models is that in the former the Coulomb interaction introduces a nontrivial 
local dynamic that is preserved in the mean-field approach. The molecular field for the electronic 
mean-field theory is not a constant number, but is instead a function of the energy, hence the 
name "dynamical mean-field theory". These local electronic mean-field theories have a rich 
history, dating back to Migdal and Eliashberg's theory of superconductivity |33| ; Metzner and 



Vollhardt |£J were the first to realize that these mean-field theories become exact solutions in 
the limit of infinite dimensions. 

Although the solution of the effective local system is still a highly nontrivial matter, it 
has the great advantage that one is able to work directly in the thermodynamic limit. Note 
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that the lattice structure has not been completely eliminated from the problem: It enters 
indirectly via the free density of states (DOS), when momentum summations are converted to 
energy integrals, and it enters in the evaluation of susceptibilities, which still maintain a (weak) 
momentum dependence. 

Two methods have proven to be most successful in solving the remaining local problem: a 
QMC scheme based on the work of Hirsch and Fye @, |35], [3(| and a perturbational method 
known as the non-crossing approximation (NCA) |36[. Both methods have their merits and 
disadvantages and may be viewed as complementary to each other, providing a means of ac- 
cessing many different quantities of interest. A more detailed discussion of the limitations of 
both approaches can be found in our earlier publications 
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The knowledge of the single-particle self energy and the local two-particle vertex functions 
(which can also be determined from the same local problem) enables one to calculate physical 
quantities. As an illustration of this, an extended discussion of magnetism in the Hubbard 
model may be found in Refs. |3^, |38], [39|. Other quantities of interest include the transport 
coefficients, the optical conductivity, the thermopower, and the Hall coefficient. For example, 
the conductivity can be calculated exactly in the dynamical mean-field theory by the following 
procedure: Figure |3] shows the leading diagrams in the expansion of the conductivity. As 
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Figure 3: Diagrammatic representation of the first two contributions to the conductivity. The 
second diagram contains a full particle-hole vertex insertion. The latter is momentum indepen- 
dent, i.e. the k-sums on both sides can be performed independently. Since the current vertex and 
the single-particle Green functions are of different parity (with respect to their k-dependence) 
the second and all higher-order diagrams identically vanish. 

was discussed earlier, the particle-hole vertex appearing in the second and all higher-order 
diagrams is momentum independent. This means that the /c-sums on the left and right end of 
all diagrams, except for the simple bubble, may be performed independently. Since the current 
operator contains the fc-gradient of the kinetic energy, and the Green functions are /c-dependent 
only through e-r, these sums identically vanish. Thus the conductivity is given by the simple 
bubble only [|17], |4(| . The evaluation of the bubble leads to 

"-<"> = Wa S *^-7 -N^{WJ M<)M< + ") , (5) 

— oo kcr x ' 
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which reduces to 



a 




(6) 



for the dc conductivity, where a is the lattice constant, the spectral weight satisfies At 



-Im 



, and /(e) := 1/[1 + exp j3e] is Fermi's function. With h/e 2 « 2.6 • 10 4 S1 the 
constants in front of Eq. (^j) can be evaluated to yield cx ~ 10~ 3 . . . 10~ 2 [(iiQcm)^ 1 }. 

The Drude weight D, may be determined by extrapolation of the Matsubara- frequency 
current-current correlation function using the method proposed by Scalapino et al.||l]|. This 
method sets a criteria to determine if the ground state of a system is a metal, insulator, or 
superconductor, by determining the asymptotic form of the current-current susceptibility in 
the x-direction, A xa; (q, iv n ) where u n = 2nirT. Specifically, the Drude weight, D, is given by 



D 



7r lim 



[-T x )e 2 -A xx (q = 0,2mT) 



(7) 



Here, (—T x ) is the average kinetic energy per site, divided by the number of lattice dimensions. 

Comparison of Eq. fl5|) with the standard expressions for transport coefficients |42| in the 
relaxation-time approximation, shows that a variety of other transport coefficients may be 
calculated if one identifies 



1 

N 



E 



(9) 



as the transport relaxation time. For example, the electronic contribution to the thermopower 
becomes |TJ2| 

S = -^rP^ , ^ « 86 W/K] 
\e\ Ln |e| 

where Ljk are the standard transport integrals 

df(e) 



de 



rHe)e k - 1 



(10) 



It is known from studies of heavy fermion systems that the relaxation-time approximation 
is sufficient to understand most zero- field properties. The Hall coefficient, however, is more 
sensitive to the approximations involved: Within the standard relaxation-time approximation 
the Hall coefficient satisfies 

-\e\ L 21 

Note that both L\\ and L 2 \ are positive, implying that the Hall-coefficient is always negative and 
that the transport is electronlike. Experiments, on the other hand, show that Rh can change 
sign in intermediate or low temperature regions, and that the Hall coefficient is usually positive 
with an anomalous temperature dependence. An explanation for this theoretical deficiency was 
proposed by the phenomenological introduction of a "skew-scattering" term |43|]. However, a 
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more refined treatment of the field-dependent conductivity can be performed The result 
is 



a 



.H 

xy 



27r 2 \e\ 3 aB 

~^h 2 



du 



1 



du / N 



E 



dk x 



2 d 2 e T , 



dk 2 



(12) 



for the Hall conductivity. Here, B is the external magnetic field which points in the z- 
direction. The constants in Eq. (O) can be rearranged according to \e\ 3 aB/h 2 = ao\e\a 2 B/h and 
\e\a 2 B/h m 10~ 5 B ~ . Inserting the values for cr , we obtain R H = cr xy /a 2 x ~ 1CT 8 . . . 10" 9 -B [m 3 / 
as the unit for the Hall coefficient. 

In the limit of large coordination number, and for a simple hypercubic lattice, the /c-sums 
in Eqs. (El) and (|P2| ) can be further simplified to yield 



de, 



N 



ktj 



- \ deN (e)K(e) 



(13) 



and 



dk 2 y 



2d 2 



deN (e)eK(c) 



(14) 



with K(et) an arbitrary function depending on k through only, and iVo(e) 



-(e/t*)S 



the noninteracting DOS. 

Examination of Eqs. |^|12|-fil| shows that the electronic conductivity is a 1/d effect and that 
the Hall effect enters to order 1/d 2 . Note, however, that Rh = &x y /& xx enters to zeroth order! 

Another important probe of the high-temperature superconductors involves nuclear mag- 
netic resonance (NMR). The spin- lattice relaxation time T\ characterizes the time needed to 
align the nuclear spins along the direction of the field. The local NMR-relaxation rate satisfies 



— = T lim Im ^^ 



(15) 



where x( a; ) is the local dynamic spin susceptibility. Similarly, the spin-echo decay rate (or 
transverse nuclear relaxation rate) T 2 g is another probe of the spin dynamics in the cuprates. 
T 2 g satisfies a more complicated relation Jl5|, (4(J 



■ 1 - 


2 


-T2G- 


= c 





Y.nk)\{k) 2 -{Y^F{k) 2 X {k)f 

k k 



(16) 



with x(k) the static momentum-dependent spin susceptibility, C an overall normalization factor, 
and F(k) the relevant form factor. The form factor involves the local site and the nearest- 
neighbor shell, and assumes the form 



F(k) 



2t 

1 + 7^ 



(17) 
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in infinite dimensions. Here, 7 is a constant of order 1. The static susceptibility turns out to 
be a function of e^/Vd =: X (in infinite dimensions), with X = corresponding to the local 
susceptibility. The integrals over the form factors can be performed, to yield 



" 1 " 


2 


c 


-T2G- 




~d 



27V (X = 0) 



d X {X) 



dX 



x=o 



;is) 



The inverse spin-echo decay rate is a 1/ \fd effect in infinite dimensions, and it is effectively 
proportional to the local susceptibility, since the derivative of the momentum- dependent sus- 
ceptibility with respect to X is an order of magnitude smaller than the susceptibility itself in 
the region of interest (see, for example Fig. 1 of ||38|| ). In fact, we neglect the derivative term 
when evaluating T 2 g P2G ^ VxP^ = 0)]- 



3 Discussion of the results 
3.1 Single-particle properties 

Theory. Before we present results for the transport properties, we will summarize some of 
the basic physics in the Hubbard model that is obtained from the dynamical mean-field theory. 
The Coulomb parameter U is fixed at U — 4. This choice may appear to be arbitrary at 
a first glance, however, the main effect of U in the strong-coupling regime is to determine 
the characteristic low-energy scale ||37||, and U = 4 is a convenient choice for numerical and 



presentational reasons. 

As the temperature is lowered, the Hubbard model in infinite-d is found to always be a 



Fermi liquid |3Q, 17], except for the region of phase space where it is magnetic |p8|, p9|. A Fermi 



liquid is defined by a self-energy that has the following structure: 

ReE(to + iO + ) = ReS(O) +u(l - Z) + O(to 2 ) , 

Im£(u; + i0 + ) = -T + O(co 2 ) , (19) 

with T oc T 2 for temperatures T < To the characteristic low-temperature scale. [The Fermi 
temperature T Q decreases to zero as half filling is approached||18||, and the Fermi-liquid-theory 



form of Eq. ( |19|) still holds for moderate temperatures, with the only change being T oc T for 
T > Tq.} The spectral weight, then assumes the form 



7r T 2 + [ujZ + e F - e^Y 



with the Fermi level defined by ep '■= n — ReS(O) and (u) denoting the (rather structure- 



less) incoherent contributions to the spectral function. The spectral function includes a delta 
function at zero temperature [^(cu) — > 5{ujZ + ep — 6%) + A 1 ^ 10 ^)] because the broadening 
T vanishes in that limit. The single-particle DOS, N(u), is defined to be the integral of the 
spectral function over all momentum, 

A^):=E^M • (21) 
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Figure 4: Single-particle DOS of the Hubbard model for different dopings 5=1 — (n) at an 
inverse temperature [3 = 43.2. At half filling (5 = 0) the DOS has a gap between the lower 
and upper Hubbard band. Away from half filling a resonance occurs at the chemical potential in 
addition to the lower and upper Hubbard bands. This resonance becomes broader with increasing 
5 and finally merges with the lower Hubbard band. The dotted curve at 5 = displays the 
(shifted) noninteracting DOS. 



Figure |] shows the evolution of the single-particle density of states as a function of doping 
for a fixed temperature (/3 = 43.2). In addition, the (shifted) DOS for the noninteracting system 
has been added for comparison (dotted line at 5 = 0). At half filling, the doping satisfies 5 = 0, 
and the system shows a pseudogap in the DOS. Note that the lower and upper Hubbard bands 
centered at uj « ±U/2 have the same approximate width as the free DOS, but are decreased 
by a factor of two from the unperturbed height due to the correlations. Away from half filling 
(5 ^ 0), a sharp resonance appears near the chemical potential. As the doping 5 increases, the 
width of this resonance also increases and it starts to merge with the lower Hubbard band. For 
dopings 5 > 0.4 both low-energy peaks become indistinguishable, implying that the system has 
become an uncorrelated metal. The upper Hubbard band seen in Fig. [| is, on the other hand, 
well separated from the low-energy excitations by a pseudogap of order U and thus contributes 
only to high-energy features in e.g. the optical conductivity JT7|. 



Figure [| collects the behavior of the DOS (Fig. |5|a) and the real and imaginary parts of the 
single-particle self energy (Fig. |5|b) as the temperature is reduced at fixed doping (5 = 0.188). 
The resonance at fi appears to be strongly temperature dependent, vanishing as T increases. 
For low temperatures, T <C T , the peak is centered near the chemical potential, and the height 
saturates to the value of the noninteracting DOS N (u = 0) = I/^/tt ~ 0.56. The upper bound 
for the local DOS follows from the following simple argument (that does not require Fermi- 
liquid-theory behavior): Since the self energy is k- independent, the interacting DOS may be 
written as 

oo 

N(u) = J deN (e) Im— (22) 

— oo 

with ( := uj — E(u; + i0 + ) being a complex number. Applying the mean- value theorem, we 
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Figure 5: (a) The evolution of the density of states when U = 4 and S = 0.188. The development 
of a sharp peak at the Fermi surface is correlated with the reduction of the screened local moment 
TxuiT), as shown in the inset. Hence the development of the peak may be associated with a 
resonant Kondo screening of the spins. 

(b) The real and imaginary parts of the self energy for various temperatures when U = 4 and 
5 = 0.188. Note that as the temperature is lowered, Im£(co>) becomes parabolic in u, indicating 
the formation of a Fermi liquid. 



obtain 

oo 

N(u) = jV (O J de (--) Im-— = jV (0 < N o {0) (23) 

— oo 

with a suitably chosen £. This implies that we must always have N(u) < N (0) = I/^/tt, 
but it does not tell us whether and where this maximal value will be reached (although the 
maximum is attained at uj = — ep/Z as T — > for a Fermi liquid). Since, on the other hand, 
the local DOS is obtained from an effective single-impurity Anderson model, we expect strong 
resonant scattering at [i in the correlated limit, i.e. the system tends to provide a large DOS 
near the Fermi level. However, it cannot build up DOS everywhere, but has to adjust itself in 
such a way that the sum rule in Eq. (|23| ) is fulfilled. The system, therefore, satisfies a delicate 
self-consistent Friedel's sum rule. 

The development of the observed resonance as T — ► is accompanied by a strong reduction 
of the effective local magnetic moment Txu{T) (see inset to Fig. ||a). This simultaneous ap- 
pearance of the resonance at fi and vanishing of the effective local moment suggests that this 
effect may be attributed to a Kondo-like quenching of the local magnetic moment in the system 
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Figure 6: (a) Angle resolved spectra for U = 4, /3 = 43.2 and (n) = 1 (a) and (n) = 0.8122 (b) 
along the Y-M direction of a two-dimensional Brillouin zone. As expected at half filling, the 
system shows a gap at fi throughout the entire zone. Away from half filling a quasi-particle peak 
develops when the Fermi surface is crossed. In addition, there are strongly damped features at 
the position of the upper Hubbard band. The data for half filling were obtained with the NCA. 



18, 37 1 . The corresponding low-energy scale, which also sets the width of the resonance at jj 



seen in Figs. f| and j|, may be estimated from l/xu(T = 0) and is tabulated in Table 1. Note, 
however, that unlike conventional Kondo systems described by Anderson's model, where one 
uncorrected band couples to a localized correlated state, leading to the quenching of the lat- 
ter's moments, there exists only one type of electron in the Hubbard model (|l|) which provides 
the band that is used to quench its own moments. This is a new effect which we propose to be 
called the "collective single-band Kondo effect" . 

We can also examine the angle-resolved photoemission spectra (ARPES) obtained from the 
self-energy. Although it is a crude approximation, we assume that the local self energy describes 
the self energy of a two-dimensional Hubbard model on a square lattice. Taking a cut along the 
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Figure 7: The data in Fig. $b for T-M multiplied by Fermi's function (at (3 = 43.2) and 
restricted to the low- energy part. This is the type of spectrum that would typically be observed 
in an ARPES experiment (see e.g. WM)- Data for k beyond 0.4(7r, 7r) were dropped because 
these features are cut off by Fermi's function. 



S = l-(n) 


T 


0.0680 


0.0177 


0.0928 


0.0273 


0.1358 


0.0478 


0.1878 


0.0730 


0.2455 


0.1074 



Table 1: Values of the Kondo scale Tq for different dopings when U = 4. 



Y—M direction in the Brillouin zone we obtain the spectra shown in Fig. |]for U = 4, (3 = 43.2, 
(n) = 1 (Fig. |5|a) and (n) = 0.8122 (Fig. ^b) (note that this approximation does not distinguish 
between a cut along the diagonal r — M or along the k x and k y axes V — X — M because the self 
energy is local, and does not depend upon k). While the system is clearly insulating and non 
Fermi-liquid like at half filling, a typical quasiparticle peak appears when the Fermi surface is 
crossed away from half filling. In addition two strongly damped bands (the lower and upper 
Hubbard bands) can be seen at higher energies in both figures. When the data of Fig. |b is 
multiplied by Fermi's function (at f3 = 43.2), one arrives at the typical ARPES result of Fig. |7|. 



This is quite analogous to what is seen in QMC simulations for true 2D-clusters [[47 



Using the peak positions in Fig. |B| as a definition of the quasiparticle energies E%, yields 
the bandstructure shown in Fig. || At half filling (Fig. |]a) there are two cosine-like bands 
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Figure 8: Bandstructure for the Hubbard model obtained from the self-energy in Fig. [6|. 



below and above /i, which represent the lower and upper Hubbard bands. They are separated 
by an indirect gap of order U/2. Off half filling, the upper and lower Hubbard bands are 
flattened relative to their values at half filling (with the most flattening occurring near the 
Fermi level). In addition, a dynamically generated flat quasiparticle band can also be seen. 
This band is the flattest of the three and has surprisingly small dispersion near the X point, 
which compares qualitatively to those obtained for small two-dimensional Hubbard-clusters 
[ fT7| , pgfl , and supports the conjecture that the dynamical mean-field theory already gives an 
accurate picture for the single-particle properties down to d — 2. Note that the flatness near 
the X point arises from both the flatness of the noninteracting bandstructure, and the many- 
body renormalizations that reduce the quasiparticle bandwidth. 
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Experiment. Only qualitative comparisons can be made between experimental photoemis- 
sion spectroscopy (PES) and inverse photoemission spectroscopy (IPES) and theoretically gen- 
erated spectra, because the experimental resolution broadening, selection rules for the matrix 
elements, and the internal dynamics of the scattering and relaxation processes, ultimately al- 
ter the experimentally measured spectrum from its theoretical counterpart. In addition (with 
the exception of first-principles calculations), simplified electronic models are used, that are 
restricted to specific energy regions (mostly the low-energy ones). Nevertheless one may at 
least identify certain features of the model with trends found in experiments. 

The measured PES and IPES spectra of the cuprates have as their most interesting trend, 
an increase of spectral weight close to the valence band of the insulating parent compound as 
the system is doped (for a review see ]|{J ) . This behavior is completely different from what one 
expects in a simple rigid-band picture in which the spectrum should be insensitive to doping 
with only the Fermi energy changing. Obviously, a similar trend is observed in the theoretical 
spectra where the quasiparticle resonance due to the collective Kondo effect develops right at 
the top of the lower Hubbard band as the system is doped (see Fig. |J). 

The quasi-two-dimensional nature of the cuprates is especially appealing to the application 
of ARPES measurements. Accordingly the cuprates have been exhaustively studied with this 
technique pD , [51], |53|, |5|J. The general behavior found for the low-energy portion of the 



ARPES results agrees with our results in that an experimental feature that might be iden- 
tified with the quasiparticle peak |2(J crosses the Fermi level in the T-M direction just like 
the theoretical result in Fig. [j]. Assuming the low-energy excitations to be Fermi-liquid-like 
[ |55|1 , produces a quasiparticle bandstructure that also shows flat bands near the X point of 
the Brillouin-zone 0, [57j identical to the theoretically generated quasiparticle band in 
Fig. H Since our results were obtained within a dynamical MFT that is rather insensitive to 
details of the underlying lattice structure, the agreement between experiment and theory is 
strong evidence that the low-energy single-particle dynamics are produced by strong electronic 
correlations independent of the dimensionality! 

3.2 Optical conductivity. 

Theory. The optical conductivity a(w) is another important probe of a strongly correlated 
system. It measures the rate at which electron-hole pairs are created by photons of frequency to. 
Figure || shows the results |17|, [58|] for a{uj) obtained from Eq. fl5|) when 5 = 0.068 for a variety 



of temperatures (a), and results for (3 = 43.2 and various dopings S (b). The Drude peak at 
uj = develops with decreasing temperature. In addition there appears a small mid-infrared 
peak at to ~ 1. As shown in Fig. ^b, this peak is more pronounced for small 5 but remains 
visible even at larger doping. In addition, it is strongly temperature dependent and clearly 
visible only for the lowest temperatures. We attribute it to excitations from the lower Hubbard 
band to the quasiparticle band at the chemical potential. The last feature in <j(uj) is a roughly 
temperature-independent peak at uj w U due to the charge excitations from the lower part of 
the spectrum, i.e. from the lower Hubbard band and the quasiparticle peak at fi, to the upper 
Hubbard band (cf. Fig. H). 

The insets to Fig. |9| show the development of the Drude weight as obtained from Eq. (|7|) 
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Figure 9: (a) Optical conductivity vs. uo for various temperatures when 5 = 0.068. Note that 
at low temperatures, when the Kondo peak becomes pronounced in the DOS, a mid-IR feature 
begins to appear in o~(uS). As shown in the inset, if we fit this data to a Drude form, then the 
width of the Drude peak is found to increase roughly linearly with T . 

(b) Filling dependence of the optical conductivity when U = 4 and j3 = 43.2. Note that for 
larger 5 the mid-IR and Drude peaks begin to merge, so that the latter is less distinct. The 
inset shows the evolution of the Drude weight D as a function of doping. D is determined by 
the extrapolation method of Ref. j^J^j. 

as a function of doping (Fig. |9|b) and the width of the Drude peak as function of temperature 
(Fig. [|a). The latter was obtained by fitting the generic form 

D r , . 

g (fa,-»0 =- - (24) 

7T 1 + T z UO A 

to the low-frequency regions in Fig. |9]a. The Drude weight initially increases linearly with 
5, then saturates to its maximal value at 5 ~ 0.5, before decreasing. This behavior can be 
understood in terms of a simple picture: The Drude weight is determined by the carrier density 
and effective mass via D oc n/m*. From the doping-dependence of the quasiparticle peak in 
the spectra in Fig. |] one may assume m* _1 oc 5. The carrier density, on the other hand, is 
given by n ~ 1 — 5, i.e. D oc 5(1 — 5). This expression explains the behavior for small doping 
as well as the maximum, which should lie at 5 max ~ 0.5, and explains how the character of the 
carriers changes from being holelike near half filling to being electronlike at low density. 

The width of the Drude peak, 1/r, displays a linear behavior 1/r ~ T for T < 0.1. This 
dependence may be traced back to the development of the Kondo peak below T . Note that 
the intercept of the linear region does not lead to 1/r — » as T — >• as required, implying that 




17 



for very low temperatures the Drude width decreases more rapidly (1/r ~ T 2 ), as expected for 
a Fermi liquid. 



Experiment. When one tries to make contact with experiment, care has to be taken about 
the energy scales. The single-band Hubbard model should only be used to describe low-energy 
features of the cuprates It does not make sense to use our results beyond uj ~ 2 in 



a comparison to experimental data, because the higher-energy bands, corresponding to the 
charge-transfer insulating behavior of the parent compounds, have been neglected in the map- 
ping to a single-band Hubbard model. 

Two types of experiments have been performed to determine the optical conductivity of 
the cuprates: A Kramers- Kronig analysis was employed to extract o{uj) from reflectivity mea- 
surements (|, m |6l|, f^, U; and photoinduced absorption has also been used |)3|, |64|. These 
experiments yield five low-energy trends: (i) The mid-IR peak maximum moves to lower fre- 
quency and merges with the Drude peak as the doping increases; its spectral weight grows very 
rapidly with doping; (ii) At a fixed value of doping, spectral weight rapidly moves to lower 
uj as T — > 0, but the total weight in the Drude plus mid-IR peaks remains approximately 
constant; the width of the Drude peak decreases linearly with T; (iii) The insulating phase 
has a charge-transfer gap; upon doping, the optical conductivity initially increases within the 
gap region; (iv) There is an isobestic point, or nearly isobestic behavior (in the sense that the 
optical conductivity is independent of doping) at a frequency that is approximately one half 
the charge-transfer gap; (v) more than one peak is observed in the mid-IR region. 

Most of these trends are observed in the theoretical model: the mid-IR peak is observed to 
move to lower frequency and merge with the Drude peak as a function of doping; the optical 
conductivity rapidly increases within the gap region as the system is doped; and there is an 
isobestic point at a frequency uj « 2. In addition, spectral weight is transferred to lower 
frequencies as the temperature is decreased, and the Drude width depends linearly on the 
temperature for a wide range of T, but the total Drude plus mid-IR spectral weight increases 
as T — > 0. The theoretical model also does not display multiple mid-IR peaks. 

These experimental features in the optical conductivity are usually attributed to either 
phonons or impurities, but judging from our results, the low-energy feature may also be due 
to excitations from the lower Hubbard band to a dynamically generated quasiparticle band at 
the chemical potential, i.e. connected to local or short-ranged spin-fluctuations. Let us stress, 
though, that we do not want to make the point that one is able to explain the normal-state 
of the cuprates in all respects quantitatively by studying the single-band Hubbard model in 
the dynamical mean-field theory. We do think, that our calculations display the underlying 
physics that drives the anomalous features found in the normal state of the cuprates, and they 
determine the order of magnitude of the corresponding temperature scale. 



3.3 Transport coefficients 

Theory The creation of a dynamical low-energy scale, like the Kondo temperature observed 
in the Hubbard model above, is known to be accompanied by interesting and anomalous features 
in transport coefficients. In our previous studies we already described the resistivity and NMR 
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relaxation rate as striking examples for this behavior. Figure [H] summarizes these results. 
Most prominent is the pronounced linear region in p(T) which increases with increased doping. 
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Figure 10: (a) Resistivity versus temperature for several different dopings when U = 4. The 
open (filled) symbols are for NCA (QMC) results. There is essentially exact agreement between 
the NCA and QMC for high temperatures. The slope in the linear regime, determined by a 
linear least squares fit, increases linearly with 1/5 as shown in the inset. The units on the 
vertical axis are 10 3 'Qcm. 

(b) NMR relaxation 1/Ti vs. temperature for different dopings at U = At* . The lines are linear 
fits in the anomalous region. 



The slope of this region is proportional to 1/5 as shown in the inset to Fig. |10|a. The thermal 
conductivity k was also calculated; however it is not plotted, since to a very good approximation, 
it follows the Wiedemann-Franz law k oc T/p. Figure |T0| b presents the results for the NMR 
relaxation rate (T\). Here, too, a rather anomalous variation with both temperature and doping 
is found. A linear region in 1/T\ develops as the doping increases, and the slope changes sign 
for 5 ~ 0.15. The doping dependence of 1/Ti is reduced as the temperature increases, but does 
not disappear at the temperatures that can be reached by the numerical analytic continuation. 

It has been recently argued that the constancy of the ratio T{T '/T 2 q is a test for the quantum 
critical region of the two-dimensional Heisenberg spins in the Cu-0 planes of the cuprates 
|65j| (even though this theory neglects the charge degrees of freedom, assuming that the most 



important effect of the holes is to add disorder into the spin system). However, as seen in 



Fig. 11, this ratio also becomes flat at intermediate temperatures within the dynamical mean- 



field theory in infinite dimensions. The constant value of the ratio T1T/T2G increases with 
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Figure 11: Ratio of the longitudinal (7\) and transverse (T 2g oc 1/xioc) NMR-relaxation rates 
multiplied byT. Note how this ratio becomes flat at intermediate values ofT, which is supposed 
to be a signal for the quantum- critical regime. The constant value of this ratio increases with 
doping, but is nearly constant near half-filling. 



doping, but does not change significantly for the lowest values of the doping. 

Finally we want to report calculations of two other interesting transport coefficients, namely 
the thermopower (S) and the Hall coefficient (Rh)- The strange behavior of the latter gave 
rise to a number of speculations about different scattering mechanisms for transport with and 
without magnetic field [^] and it is thus interesting to look at the behavior of this quantity 
within our scheme. Fig. |l| compiles the results for S (Fig. |T2]a) and Rh (Fig. 0b) for four 
different dopings as function of temperature. The thermopower shows the sign change at 
intermediate temperatures that is characteristic of correlated materials. The Hall coefficient, 



in Fig. |12|b, is more interesting. It is positive for high temperatures, and displays a maximum 
at intermediate temperatures, followed by a strong decrease for lower T, eventually becoming 
negative (qualitatively similar features have been seen in simulations of the two-dimensional 
Hubbard model Interestingly, the position of the maximum is weakly sensitive to doping 

while, as shown in the inset, its height roughly decreases with increasing doping like Rh ~ 1/5. 
Finally, we show the quadratic behavior of the Hall angle (cot 0# = p xx / p xy = l//i, where /z is 
the charge carrier's mobility). The Hall angle is calculated with the NCA so that a dense set of 
points may be presented at high temperatures, where the analytic continuation of QMC data 
becomes numerically expensive. Here the agreement between the NCA and the QMC data is 
essentially exact (cf. Fig. [TJ]). Fig. plots the inverse mobility as function of temperature for 
two dopings (cf. Fig. |T] for typical experimental temperature and doping dependence in the 
cuprates). The Hall angle is plotted as a function of T 2 in the inset to Fig. [13]. It shows a 
roughly linear behavior on this scale. 



Experiment. As we have pointed out the theoretical behavior is consistent with the 
normal-state properties of hole-doped high-T c compounds. As is well known, the experimental 
trends observed include Q (cf. also Fig. [I]): (i) A (sub-) linear variation of the resistivity 
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Figure 12: Thermopower (a) and Hall coefficient (b) for four different dopings 5 as function of 
temperature. In (b), the open symbols (dashed lines) are QMC (NCA) results and in (a) the solid 
lines come from a fit to the QMC data (open symbols). In the inset to (a), the Hall coefficient 
at a fixed temperature Rh{T = 0.1852) is plotted versus 1/5, indicating that Rh(T = 0.1852) 
increases roughly in proportion to 1/5, consistent with experimental results for the cupratesfiUBj]. 
The units on the vertical axis are 86/iV/K for the thermopower (a) and 10~ 9 m 3 /C for (b). 

with temperature, the slope or absolute value of p(T) for a given T decreasing with 1/5; (ii) 
The NMR-rate drops with increasing doping due to a decrease of the (local) spin fluctuations 
and shows a linear tail with positive slope at intermediate temperatures; There is also a ten- 
dency towards a change in sign of the slope of this linear region with decreasing doping; (iii) 
The spin-lattice relaxation rate T\ becomes doping independent at intermediate values of the 
temperature, and the ratio T{T /T 2 g becomes constant over a similar temperature range; (iv) 
The Hall coefficient goes through a maximum whose position is roughly independent of doping 
but whose height decreases oc 1/5; (v) The Hall angle cot 6# is usually found to vary like T 2 
over a considerable temperature range. Since all these features (except the constancy of 1/T\ 
with respect to doping at intermediate T) are consistently obtained from the dynamical MFT, 
including the peculiar behavior of the Hall angle (without any necessity to resort to an exotic 
ground-state or scattering mechanism) and the constancy of the ratio T1T/T2G, we anticipate 
that these anomalies are intrinsic properties of strong local correlations and do not specifically 
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Figure 13: Inverse mobility \ j\L — p/Rh calculated with the NCA for two dopings as function 
of T for large T. The NCA was used since it is essentially exact at these temperatures (cf. 
Fig. and significantly less computationally expensive, thus allowing us to present a dense 
data set. In the inset the Hall angle cotO# is plotted as function of T 2 . Note the apparent 
linear behavior. 



depend on the low dimensionality of the cuprates. This statement is further supported by the 
observation that similar features in the NMR and EPR relaxation rates and the Hall coefficient 
are also found in "conventional" heavy-fermion materials [16, 69 . 



4 Summary and conclusion 

We presented results for a variety of transport properties that support the conjecture that most 
of the anomalous normal-state properties of the high-T c compounds are intrinsic properties 
resulting from correlation effects in a single-band model. These anomalies include a linear 
in T resistivity and anomalies in the temperature dependence of the NMR-relaxation rate 
and thermopower. We also presented a method to calculate the Hall-conductivity and Hall 
coefficient that does not rely on the (conventionally used) relaxation-time approximation or 
the introduction of a new scattering mechanism like Coleman's "skew-scattering". This new 
method enabled us to obtain physically sensible results for the Hall coefficient, showing a 
positive sign and a temperature dependence characteristic of strongly correlated systems. We 
also found that the Hall angle or inverse mobility of the carriers shows a clear T 2 -behavior 
as also observed in the cuprates. The optical conductivity was found to have a Drude peak 
at small u>, a charge-excitation peak at uj ~ U, and a mid-IR feature that is attributed to 
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excitations from the lower Hubbard band to the quasiparticle band dynamically generated at 
the Fermi level. 

Most of these properties also show a distinctive dependence upon doping: The Drude weight 
increases linearly with doping 5 and both the slope of the linear region in the resistivity and 
the Hall coefficient increase ~ 1/5. The mid-IR feature in the optical conductivity initially 
increases rapidly with doping, and then decreases as 5 increases further (but is still visible up 
to 5 « 0.20). 

These features are also found in the high-T c materials and to some extent in heavy-fermion 
or mixed-valence compounds. From our results it seems that the appearance and overall tem- 
perature dependence of these anomalies is due to the existence of strong local correlations 
that lead to a dynamically generated (strongly temperature-dependent) low-energy scale aris- 
ing from a Kondo-like screening of the (local) magnetic moments. The anomalous regions are 
more pronounced in the present case, which may be attributed to the fact that in the Hubbard 
model only one band exists, i.e. the electrons that form the moments are also responsible for 
their screening, while these different tasks are split between at least two bands in the periodic 
Anderson model. To distinguish between these different physical situations we suggest that the 
singlet formation found in the Hubbard model should be labeled as a collective "single-band 
Kondo-effect" . 

Motivated by these considerations, we propose that the peculiar anomalies found in the 
cuprates may be viewed as a crossover phenomenon from a high-temperature "normal" phase 
to a renormalized Fermi liquid as T — > 0. The anomalous temperature dependence of physical 
quantities is then obtained from the peculiar T-behavior of the developing quasiparticle band 
at [i. Note that in several heavy-fermion compounds true Fermi-liquid behavior is not observed, 
because the development of the Fermi liquid is preempted by a phase transition into an ordered 
state. It was fortuitous that the first heavy-fermion compounds studied did have transition 
temperatures much smaller than Tk, enabling one to observe the formation of the heavy Fermi 
liquid first. In the cuprates, on the other hand, it may be that the relevant Kondo scale is of 
the same order as the transition temperatures, and we are in a situation where the mentioned 
crossover is observed, but that the system undergoes a phase transition before a Fermi liquid 



forms W0\ 
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